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It is shown that the Hilbert-Schmidt (HS) norm and distance, unlike the trace norm and distance, 
are generally not contractive for open quantum systems under Lindblad dynamics. Necessary and 
sufficient conditions for contractivity of the HS norm and distance are given, and explicit criteria in 
terms of the Lindblad operators are derived. It is also shown that the requirements for contractivity 
of the HS distance are strictly weaker than those for the HS norm, although simulations suggest that 
non-contractivity is the typical case, i.e., that systems for which the HS distance between quantum 
states is monotonically decreasing are exceptional for N > 2, in contrast to the case N = 2 where 
it is always monotonically decreasing. 



I. INTRODUCTION 

The trace norm and induced trace distance play an im- 
portant role in quantum information theory. One of its 
most important features, widely used in research papers 
and popular textbooks [TJ[2], is its contractivity for trace- 
preserving quantum evolution, first proved by Ruskai |3J. 
It was once conjectured that contractivity would extend 
to other more intuitive norms such that the Hilbert- 
Schmidt (HS) norm |U [5], but flaws in the original ar- 
gument were soon discovered and an explicit counter- 
example of a trace-preserving map with non-contractive 
HS norm was provided by Ozawa [B]. The more general 
question of when a positive trace-preserving [PTP] map 
between matrix spaces is contractive with respect to the 
p-norm for p > 1 was recently considered in [7] , where it 
was shown that PTP maps are contractive for p > 1 in 
general if only if they are unital. This does not answer 
the question, however, when a PTP map defined on a 
subset of a matrix space is contractive with respect to 
a particular p-norm such as the HS norm, and contrac- 
tivity depends on the subspace. In particular, the HS 
norm restricted to the trace-zero hyperplane in the space 
of Hermitian matrices may be contractive while it is not 
on the whole matrix space. 

Contractivity of the distance between quantum states 
in the sense that the distance between two quantum 
states is monotonically decreasing, is of particular inter- 
est in the area of open system dynamics. We consider 
under what conditions the evolution of an open quan- 
tum system subject to semigroup dynamics governed by 
a Lindblad master equation is contractive with respect 
to the HS distance. It can be shown that for N — 2 both 
the HS norm and distance are always contractive, while 
for N > 2 both are contractive only for small subsets of 
open systems, with the set of open systems for which the 
HS norm is contractive being strictly smaller than the set 
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of open systems for which the HS distance is contractive. 
Aside from the relevance to quantum information, e.g., in 
the construction of entanglement measures, one impor- 
tant implication of the non-contractivity of the HS norm 
is in quantum control, where the HS distance between 
quantum trajectories is a common choice of a Lyapunov 
function, e.g., for local optimal control [8j [9j [10J [11] . As 
monotonicity is a prerequisite for a Lyapunov function, 
it means that the HS distance is generally not a suitable 
candidate for a Lyapunov function for open systems. 



II. HILBERT SPACE NORMS AND 
CONTRACTIVITY 

Let TC be a Hilbert space with dim H = N and let 93 [H] 
be the bounded operators on TC, and D[H\ be the trace-1 
positive operators on TL. It is not difficult to check that 

||A|| TR =iTr\/IU (1) 

is well-defined for any A £ 93 [7i] as A^A is positive, and 
satisfies the axioms of a norm, and induces a metric, the 
trace distance, in the usual way cZtr(A B) = \\ A — S|| X r. 

Among the nice features of the trace norm is contrac- 
tivity under trace-preserving maps [3], i.e., given two den- 
sity operators P\,P2 £ D[H] and a map £ : 93 [H] — » 93 [H] 
with Tr(£(A)) = Ti(A), then 

d TR (£(pi),£(p2)) < d TR ( Pl , P2 ). (2) 

This implies in particular that any (super) operator that 
maps density matrices to density matrices cannot in- 
crease the trace distance between two states. Further- 
more, if £t is a quantum dynamical semi-group then 
<^TR,(^t(pi)j £t(P2)) is monotonically decreasing, although 
monotonicity is not always strict. Indeed, in the special 
case of unitary evolution the trace distance remains con- 
stant, d.TR(£t(pi),£t(p2)) = dTR(pi,P2) for all t. 

One drawback of the trace distance, however, is that it 
is not the most intuitive distance measure for quantum 
states. A more natural choice is the Hilbert Schmidt 
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(HS) distance d KS (A,B) = ||A - 5|| H s where ||A|| H s = 
y/(A\A) is the norm induced by the HS inner product 

(A\B) m = Tr(AtS). (3) 

For Hermitian operators A = A\ B = the HS inner 
product and distance simplify to (A|j4)hs — Tr(A 2 ) and 



d HS (A,B) = ^/Tr[(A-Bn 



(4) 



If we choose an orthonormal basis {<Jk}k=i f° r the Her- 
mitian operators on Tt then the coordinate vector a — 
( a k)k=i w hh a,k = Tr(ofcA) for any Hermitian operator A 
on 7i is a vector in R w , and noting that Tr(<7fc<7^) = 
by orthonormality, shows that the HS inner product re- 
duces to the standard Euclidean inner product in M. N 



N 



N- 



(A\A) 



HS 



^2 a k a e Tr(cr fe cr £ ) = ^2 a 



k,e=i 



fc=i 



If we choose an orthonormal basis with the last basis vec- 
tor a N 2 = ^t=I, where I is the identity on Ti., then all 
Hermitian operators with constant trace Tr(v4) = c cor- 
respond to hyperplanes in M. N with a^2 = ^=c. There- 
fore, if we are only interested in Hermitian operators with 
constant trace, such as the class of density operators on 
H, it suffices to consider the reduced coordinate vector 
a = (oi, . . . , Ojv=>-i) and 

(5) 



(A\A) w =N- 1 Tc(A) + {a\a). 



Furthermore, if A, B are two Hermitian operators with 
Tr(A) = Tr(B) then Tr(A - B) = and thus 



d HS (A, B) = y/(a-b\a-b) = \\a - b\\ 



(6) 



i.e., the HS distance between two Hermitian operators of 
the same trace class, is simply the Euclidean distance of 
their associated reduced (real) coordinate vectors, which 
generalize the Bloch vector for N — 2. This equivalence 
of HS distance between density operators and the Eu- 
clidean distance between their reduced coordinate vec- 
tors makes the HS distance a very intuitive and useful 
distance measure. 



III. CRITERIA FOR MONOTONICITY OF HS 
NORM AND DISTANCE 



For N = 2, i.e., a single qubit, it is easy to show that 
the trace and HS distance agree up to a constant factor, 
and hence contractivity of one implies contractivity of 
the other. Let p\ and p2 be two qubit density operators 
and let r\ and r 2 be their respective coordinate vectors 
in K 3 with respect to the (orthonormal) basis {o^\ =l 



defined above. Then pj = J2k=i r jk^k 
and setting <x = (oi, cr 2 , 03) 



for j = 1,2 



k=l 



(7) 



Choosing er to be, e.g., the standard (normalized) Pauli 
basis, it is easy to verify that the eigenvalues of p\ — P2 are 
i^H?"! — r 2 ||, and noting that for Hermitian matrices 

Tr I A\ is the sum of the absolute values of the eigenvalues 
of A, we have 113 



<hB.(pi,P2) = \ Tr \Pi — Pal 

1 2 11 1 
= 2^lln-r 2 | 



V2 



dus(pi,p2)- 



(8) 



Hence, the HS distance for a single qubit is monotonically 
decreasing under completely positive maps and quantum 
semigroup dynamics. This argument does not generalize, 
however, as the relationship between the eigenvalues of a 
positive map and the Euclidean norm of the coordinate 
vector that was used is very specific to N = 2. 

For N > 2 we still have ||A|| H s < II^IItr from basic 
functional analysis |12j . Thus the trace norm provides an 
upper bound on the Hilbert-Schmidt norm. Applied to 
open quantum systems subject to semi-group dynamics 
governed by a Lindblad master equation 



p(t) = -t[H,p(t)}+£ DP (t), 
with C D p(t) = E d V[V d ]p(t), V d G 35 [W] and 



(9) 



v[v d ] P (t) - v dP (t)vj - \{vlv dP {t) + P (t)vjv d ), (10) 

this means that if the trace distance between any two 
quantum states goes to zero for t — > 00, for instance, 
then the HS distance must go to zero as well, but the 
convergence need not be monotonic. 

Expanding the density operator p and super-operators 
Ch and Cd with respect to an orthonormal basis for 
the (trace-zero) Hermitian matrices on Ti, it is easy to 
show that the master equation (|9| becomes an affine- 
linear equation for the coordinate vector r G R w _1 



r(t) =Ar(t) + c, 



(11) 



where A is a (TV 2 — 1) x (N 2 — 1) real matrix and 
can be computed from the Hamiltonian and 



c G 



Lindblad generators (See Appendix B). 

Proposition 1. A necessary and sufficient condition for 
the HS norm \\p(t)\\2 of any quantum state p(t) to be 
monotonically decreasing under the open system dynam- 
ics (11) is that the Bloch equation is linear, i.e., c= 0. 



Proof. If T is a PTP map acting on the Hermitian oper- 
ators on a finite-dimensional Hubert space then ||Tp|| 2 < 
\\pW2 if and only if T is unital [7]. For our open system 
dynamics, the evolution is unital, i.e., preserves the iden- 
tity I, if and only if the Bloch equation is homogeneous, 
i.e., c = 0, as the identity I is mapped to r = in the 
trace-one hyperplane the density operators live in. □ 

Proposition 2. A necessary and sufficient condition for 
the HS distance \\p\{i) — p2(t)\\ between quantum states to 
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be monotonically decreasing under the dynamics (11) is 
that the symmetric part A + A T of the evolution operator 
A be negative definite, i.e., have no positive eigenvalues. 

Proof. Denote the set of all physical coordinate vectors 
as £>r[W], corresponding to S)[7<]. We can easily see that 
the eigenvalues of A must all have non-positive real parts 
since the dynamical system is invariant on the compact 
set 2)r[H]. Moreover, if pi(0) = p\ and /^(O) = pi are 
two initial states with coordinate vectors r\ and r2, re- 
spectively, then A(t) = r\(t) — r^t) satisfies the linear 
equation 



A(t) = A A(t) 



A(t) 



A(0) (12) 



and thus we have 
d 



dt 



d 2 m ( ri (t),r 2 (t)) = (A(t)|A(t)) + (A(t)|A(t)> 

= (AA(t)) T A(t) + A(t)AA(t) 
= A(tf(A T + A)A(t). 



As A + A T is real symmetric, its eigenvalues are real, and 
the HS distance will be monotonically decreasing, if and 
only if A + A T has only non-positive eigenvalues. □ 

Note that although we know that the eigenvalues of 
A have non-positive real parts, this does not imply that 
A + A T has non-positive (real) eigenvalues in general. 

Example 1. Consider the system p(t) = T>[V]p with the 
simple Lindblad generator 




Choosing {o~k : k = 1, . . . , 8} to be the standard orthonor- 
mal basis for the trace-zero Hermitian matrices (Al) for 
N = 3, and = vie obtain 



A = 



/-I 








1 


|V3 





\ 1\ 






1 

2 

-1 




-1 










1 

I 

2 






-1 








1 

2 


-1V3 





| 














-2 





|a/3 V3 





1 

2 


1 

2 








3 
2 





1 

2 








-1 


-|V3 





-1 2" 


V i 








1 

2 


-|V3 





2 2/ 



as well as 



^(1,0,0,1,^3,0,0, 



-1) T . It is easy to 



check that A is invertible, its eigenvalues have negative 
real parts, and the system has a unique steady state r ss — 
— A -1 c corresponding to 

2-10 

r I -1 2 -1 

-1 1 



0.114 



0.08 



» 0.06 




0.04 



0.02 



FIG. 1: Non-monotonic convergence to steady state. 



However, [A, A T ] ^ 0, i.e., A is not normal and A + A T 
has a positive eigenvalue 7 = 0.1914 with eigenvector v. 
Hence, setting A(0) = av with a > chosen such that 
r(0) = r ss + av £ £>r[7Y] gives a trajectory for which 



dt 



d 2 (0) = a 2 v T (A T + A)v = 7a 2 |H| 2 > 0, 



(13) 



and thus the distance from the steady state increases at 
least initially. Indeed Fig. \l\ shows that for A(0) = 
and s(0) £ Sk[W] the HS distance dns(i) = l|A(t)|| even- 
tually converges to — as it must as A has no eigenvalues 
with real part — but it increases initially. 



Whenever A + A T has a positive eigenvalue 7 with 
eigenvalue v then the distance between any two initial 
states pi(0) and £2(0), whose corresponding real coordi- 
nate vectors satisfy r% — r 2 = av, will increase at least 
initially. Furthermore, if the system has a steady state 
p ss in the interior of D[H] then there arc initial states 
that do not converge to the steady state monotonically 
with respect to the HS distance. This is easy to see since 
for a point in the interior, it is always possible to choose 
a > such that s(0) = s ss + av 6 £)R[7i], where S)r[7Y] 
is the subset of K. _1 corresponding to physical states, 
for any v. 



IV. SUFFICIENT CONDITIONS ON LINDBLAD 
OPERATORS FOR MONOTONICITY 

A necessary and sufficient condition for the HS norm of 
a quantum state to be monotonically decreasing is that 
C D (1) = 0. Inserting this into the LME ^ leads to 
2^,2 [Vrf, Vj] = 0, which gives the following explicit result: 

Proposition 3. The HS norm of any quantum state is 
monotonically decreasing under the open system dynam- 
ics (11) if and only if the Lindblad operators V c i satisfy 



E d [v d ,v b 



0. 
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We know that contractivity of the HS norm is a suffi- 
cient condition for contractivity of the HS distance, but 
we can derive other sufficient conditions. 

Proposition 4. The distance between any two quantum 
states is monotonically decreasing under the open system 
dynamics (11) if A is normal. 



Proof. If A is normal, i.e., [A, A T ] = 0, then there exists 
a unitary transformation U £ §>U(N 2 — 1) such that A = 
UDU f , where D is a diagonal matrix. Hence, noting 
that A is real, A T = A f = UD t U t , and thus A + A T = 



U(D + D^)U^ shows that the eigenvalues of A 
twice the real parts of those of A, i.e., non-positive 



A T are 
□ 



For simple Lindblad equations we can further show 
that normality of the Lindblad operators is a sufficient 
condition for normality of A- 

Proposition 5. For a system governed by the purely dis- 
sipative LME p{t) = 2?[V]/j(i) the superoperator A + A T 
is normal if V is normal. 

Proof. Let a^ k be the (£, fe)th component of Ay- From 
the definition of the Bloch equation it follows 



= TtO^atVaj) - \ Tt(0V{o*, (14) 

where A,B = AB + BA is the anticommutator, and thus 
the (£, k)th component of Ay Ay is equal to ^ . agja k j 

and the (£, fc)th component of Ay Ay is Ylj a jt a jk- If 
Hi and Hi are two Hermitian matrices and at least one 
has zero trace, we have the identity 

MHio- 3 ) TriHzo-j) = Tv(HiH 2 ), 



and thus normality of Ay is equivalent to 

= Tr(a £ Va k V^[V, V*]) + Ti{a h Va t V^[V, V*\) 
+ Tr(Va e V^a k [V, V*]) + Tv{Va k V^ a t [V, V*[) 
+ Tr(VV*o- e VV*a k ) - Tr{V^ VaiV^Va k ), 



for all k,£, which is satisfied if [V, V'] = 0. 



□ 



erator A + A is negative definite if all Lindblad opera- 
tors Vd are normal. This is of course consistent with the 
previous observation that 5Z d [Vd, V^f] = implies con- 
tractivity of the HS norm, and hence contractivity of the 




FIG. 2: Three-level A-system with (independent) decay from 
the excited state to either of the two ground states. 

HS distance. We might thus conjecture that the two suf- 
ficient conditions for contractivity of the HS distance are 
equivalent, but this is not the case. 

Example 2. Consider the a three-level system with H = 
and two Lindblad operators 



<0 V 
Vi = ( 
,0 0, 



'0 (F 
V 2 = I 1 

,0 1 0; 



(15) 



If there are multiple Hamiltonian and Lindblad terms 
then by linearity of the master equation, the superop- 
erator splits, i.e., A = Ah + J2dAv d , wnere Ah is as- 
sociated with the Hamiltonian dynamics and Ay d corre- 
sponds to the decoherence operator Vd- It is easy to see 
that Ah is real-antisymmetric, and thus Ah + A H = 0. 
Since the sum of negative semi-definite matrices is neg- 
ative semi-definite, A + A T is negative semi-definite if norm is not. 
Av d + Ay d is negative semi-definite for all d, and the lat- 
ter is the case if Ay d is normal. Thus we can conclude 
that for a system governed by the LME (|9| the superop- 



It is easy to check that [Vd,Vj] ^ for d — 1,2 but 

Ed=i iWdi ^d] = an d thus c = 0. Thus the evolution is 
unital and both the HS norm and distance are contractive, 
but the superoperator A is not normal. 

Hence, the evolution of the system may be unital even 
if A is not normal. Similarly, the superoperator A of 
an open system may be normal even if the evolution is 
not unital, showing that the condition X)dlX<i!Kl] = 
is not a necessary condition for the HS distance between 
quantum states to be monotonically decreasing. 

Example 3. Consider a three-level A system with decay 
from the excited state to the ground state as pictured in 
Fig. [J| Neglecting the Hamiltonian part, the evolution 
equation given by two spontaneous emission processes, 
characterized by Lindblad operators 



Vi 



Lt easy to check that [Vi,!^] = diag(l,0,— 1) and 
[V 2 ,Vj] = diag(0, 1,-1) and thus J2dl v d, Vj] + 0. Thus 
the evolution is not unital, and neither of the correspond- 
ing superoperators Ay and Avi are normal, but if both 
decay processes are equally likely, the off-diagonal terms 
in the sum Avi + Avi cancel, and the resulting super- 
operator A is diagonal, hence normal. Thus the HS 
distance between quantum states under this semi-group 
dynamics is monotonically decreasing although the HS 




V. MONOTONICITY OF HS DISTANCE FOR 
NON-NORMAL A, NON-UNITAL EVOLUTION 

We have shown that contractivity of the HS norm is 
a sufficient but not a necessary requirement for contrac- 
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non-contractive 


0% 


63.7% 


95.4% 


100% 


100% 


100% 


100% 


max. no. 7 > 
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11 



TABLE I: Percentage of systems governed by the Lindblad 
master equation p(t) — T>[V]p(t) with randomly generated 
V, for which the HS distance is not monotonically decreasing 
and maximum number of positive eigenvalues 7 of A + A T 
as a function of the system dimension shows a rapid increase 
in the frequency of non-contractive dynamics for increasing 
N. For each N 1000 systems were simulated. All randomly 
generated V (and hence A) were non-normal, as expected as 
normal matrices form a measure-zero set. 
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FIG. 3: Bloch superoperator A not normal but negative 
definite. The evolution is not unital and ||p(t)||2 for po = 
diag(0, 0, 1) shows that the HS norm is not contractive. The 
the distance of p(t) from the steady state p ss , however, is 
monotonically decreasing, as is the distance between any two 
quantum states. 



tivity of the HS distance, and that normality of the dis- 
sipative part of A |14j is an alternative sufficient condi- 
tion for monotonicity of the HS distance under Lindblad 
semi- group dynamics. Finally, we show that even both 
sufficient conditions together are not necessary, i.e., there 
are systems for which A is not normal and c / 0, but 
the HS distance is still monotontically decreasing. 

Example 4. Consider the purely dissipative semi-group 
dynamics given by p(t) — T>[V]p(t) with 

/0 1 0> 
V = 1 

\0 Oy 

V is clearly not normal and neither is the superopera- 
tor A, which, with respect to the standard basis (Al) for 
the trace-zero Hermitian matrices for N = 3, takes the 
explicit form 
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Furthermore, we have c / 0, i.e., the evolution is not 
unital. Nonetheless, we can check that the eigenvalues 
of A + A T are -2 ± fv^, -§ ± and -1, all of 

which are negative. (There are only five distinct eigen- 
values as the last three occur with multiplicity 2.) Thus 
the HS distance between any two states is monotonically 



decreasing, and we can easily verify that the system has a 
unique steady state p ss = at the boundary ofD[H]. 

The HS norm, of course, is not contractive. E.g., if we 
start in the state po — diag(0, 0, 1) then the HS norm 
\\p(t)\\2 = \/Tr(p(i) 2 ) first decreases and then increases 
as shown in Fig. [#[ 

However, examples when A is not normal, the evolu- 
tion not unital and the HS norm non-contractive, but 
the HS distance is still monotonically decreasing appear 
to be increasingly hard to find in higher dimensions. Nu- 
merical tests with randomly generated simple Lindblad 
generators V suggest that when the dissipative part of 
A is not normal then there is a high probability that the 
symmetric part A + A T will have at least one positive 
eigenvalue. Moreover, both the probability of a posi- 
tive eigenvalue and the number of positive eigenvalues of 
A + A T appear to increase with the system dimension 
(See Table I). 



VI. CONCLUSION 

Unlike the trace distance, and perhaps contrary to in- 
tuition, the HS norm is generally not contractive under 
positive trace-preserving maps, except for N = 2. For 
open systems governed by a Lindblad master equation 
the necessary and sufficient conditions for contractivity 
of the HS norm translate into a necessary and sufficient 
condition for the Lindblad generators. This condition is 
also sufficient to ensure that the HS distance between 
quantum states is monotonically decreasing, but it is 
not necessary. We derive alternative necessary and suf- 
ficient conditions for monotonicity of the HS distance 
in terms of the spectrum of the symmetric part of the 
super-operator, and show that they lead to alternative 
sufficient conditions for monotonicity of the HS distance, 
which are strictly weaker than those for monotonicity of 
the HS norm. This means that the HS distance between 
any two quantum states under Lindblad dynamics can 
be monotonically decreasing even if the HS norm ||/9(i)||2 
of quantum states is not monotonic under this evolution. 
Although the criteria for monotonicity of the HS distance 
are weaker, in general it is not monotonically decreasing 
if the Hilbert space dimension is greater than 2, even 
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for systems that have a unique steady state. The non- 
monotonicity of the HS distance has important implica- 
tions for, e.g., quantum control, showing that unlike for 
Hamiltonian systems, it is not a suitable Lyapunov func- 
tion for generic open systems. It is a suitable candidate 
for a Lyapunov function only in very special cases, e.g., 
when all the Lindblad decoherence operators are normal. 
Although such systems are only a small subset of pos- 
sible systems in higher dimensions, they do include the 
important special case where the decoherence is induced 
by measurement of a Hermitian observable. 
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for 1 < r < N - 1 and r < s < N. 



APPENDIX B: BLOCH REPRESENTATION 



Let {<Jk\k=i b e a basis for the Hermitian matri- 
ces. With respect to this basis the master equation ^ 
with dissipation term (10) can be written in coordi- 



nate form as a linear matrix differential equation (DE) 
r = (L + EdD (d) ) r > where r = (r„) S R^ 2 with 
r n = Tr(pa n ) and L and are N 2 xN 2 (real) matrices 
with entries 



L mn = Tr(iH[a m ,a n ]) (Bla) 
£>W = Tr(V}a m V d a n ) - - Tr(vjV d {a m ,a n }) (Bib) 



APPENDIX A: STANDARD BASIS FOR 
HERMITIAN MATRICES 

A standard orthonormal basis for the trace-zero Her- 
mitian matrices for any N is given by {cfcl^L] -1 where 
a k = Cfe(r,s) with k — r + (s — 1)N and 

°V. = ^(|r)<s| + |*><r|) (Ala) 
tr„ = ij s (-\r){8\ + \a){r\) (Alb) 
°Vr = 7=f^ (EU l fc X fc l - r\r + l)(r + 1|) (Ale) 



where {^4, B} = AB + BA is the usual anticommutator. 

If we choose the basis such that o w2 = -4^1 and the 

v iv 

remaining basis elements form a basis for the trace-zero 
Hermitian matrices, then, noting that rjv — ^= Tr(,o) = 

is constant and thus = 0, we can define the re- 
duced so-called Bloch vector s — (t*i, . . . ,r N 2_ 1 ) T , and 
rewrite the linear matrix DE for r as affinc- linear matrix 
DE s(t) = As(t) + c for s, where A is an (N 2 - 1) x 

(TV 2 - 1) real matrix with A mn = L mn + J^d^"™- and 



//v 



i EjD W 



fN 



mN 
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